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> ■ Abstract 
O 

. We study the renormalized energy-momentum tensor (EMT) of the inflaton 

fluctuations in rigid space-times during the slow-rollover regime for chaotic 
\ inflation with a mass term. We use dimensional regularization with adiabatic 

' subtraction and introduce a novel analytic approximation for the inflaton fluc- 

■ tuations which is valid during the slow-rollover regime. Using this approxima- 

tion we find a scale invariant spectrum for the inflaton fluctuations in a rigid 
space-time, and we confirm this result by numerical methods. The resulting 
renormalized EMT is covariantly conserved and agrees with the Allen- Folacci 
result in the de Sitter limit, when the expansion is exactly linearly exponential 
in time. We analytically show that the EMT tensor of the inflaton fluctua- 
i i tions grows initially in time, but saturates to the value H 2 Hq, where H is 

(3JT), the Hubble parameter and Hq is its value when inflation has started. This 



result also implies that the quantum production of light scalar fields (with 
mass smaller or equal to the inflaton mass) in this model of chaotic infla- 



■ tion depends on the duration of inflation and is larger than the usual result 

extrapolated from the de Sitter result. 



I. INTRODUCTION 

Particle production in expanding universe, pionereed by L. Parker 0, is an essential 
ingredient of inflationary cosmology 0. The nearly scale invariant spectrum of density 
perturbations predicted by inflationary models |§ is at present inextricably related to the 
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concept of amplification of vacuum fluctuations by the geometry. The scale invariant spec- 
trum for massless minimally coupled fields during a de Sitter era was indeed computed before 
inflation was suggested. 

The calculation of the energy carried by these amplified fluctuations is then the natural 
question. To answer this question a renormalization scheme is necessary, as in ordinary 
Minkowski space-time. Ultraviolet divergences due to fluctuations on arbitrary short scales 
are common in field theory. In Minkowski space-time, infinities in a free theory are removed 
by the subtraction of the vacuum expectation value of the energy, also called normal ordering, 
the physical justification being that these vacuum contributions are unobservable. 

A similar prescription is used in order to regularize infinities in cosmological space- 
times. However, one additional problem is the absence of an unambiguous choice of vacuum, 
because of the absence of a class of privileged observers, which are the inertial observers in 
the Minkowski space-times. The idea is then to subtract the energy associated with a 
vacuum for which the effects of particle production by the time- dependence of the metric 
are minimized. This vacuum is determined by the assumption of an adiabatic expansion of 
the metric. This procedure is therefore called adiabatic subtraction In this way the 

infinities of field theory are removed. 

Even with a minimal prescription such as the above, there are several surprising effects 
accompanying the renormalized energy-momentum tensor (henceforth EMT) of a test field 
in cosmological space-times. To name a few, there could be the avoidance of singularities due 
to quantum effects 0, conformal anomalies which break classical conformal symmetries at 
the quantum level |7|], violation of the various energy conditions ||. One of the first models of 
inflation proposed by Starobinsky 15) was indeed based on the role of the conformal anomaly, 
which both avoids the singularity and produces an inflationary phase. 

While the adiabatic vacuum for a test field - and its associated EMT - can be computed 
for generic cosmological space-times, the unrenormalized EMT can be calculated analytically 
only if exact analytic solutions for the field Fourier modes are available. Because of this, 
space-times such as de Sitter have been throughly investigated |10|]5|1, since analytic solutions 
for a scalar field with generic mass and coupling are available. In the absence of analytic 
solutions, numerical schemes are implemented In a de Sitter space-time, the back- 

reaction of a test field seems important only if m 2 + £i? = 0, with m and £ separately 
different from zero, or m 2 + £i2 < [m is the mass of the test field and £ its coupling to the 
curvature R) |12|,|13J. In the former case, the EMT of the test field grows linearly in time, 



while in the latter case it grows exponentially. The EMT of a massless minimally coupled 
field is constant in de Sitter space-time |14j . 

Inflationary models based on the use of scalar fields, have an accelerated stage, usually 
called slow rollover, in which the Hubble parameter is almost frozen. During this stage it 
is rare to have exact solutions for the fluctuations. However, these inflationary models are 
more attractive than the de Sitter space-time one in furnishing phenomenological models 
whose predictions can be tested against observations. Chaotic models [jnj are the simplest 
among these. In this paper we study analytically the problem of back-reaction of inflaton 
fluctuations during the regime of slow-rollover for the case of a massive inflaton. We consider 
inflaton fluctuations in rigid space-times, i.e. we neglect metric perturbations coupled to 
them, as pioneeringly investigated by Abramo, Brandenberger and Mukhanov ||16|| . We plan 
to come back to this issue in a future work [|TT|. 
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The plan of the paper is as follows: in section II we describe the background classical 
dynamics for a massive inflaton and the novel analytical approximation for its fluctuations 
during the slow-rollover regime. In section III we discuss the normalization of quantum 
fluctuations and in section IV we compare the numerical evaluation of the spectrum with 
the analytic approximation. We discuss the EMT of inflaton fluctuations, the adiabatic 
subtraction and the renormalization in sections V-VII, respectively. In section VIII the 
problem of the back-reaction of the EMT of inflaton fluctuations is addressed and in section 
IX our novel analytic approximation is compared with the slow-rollover technique ||18|| . In 
section X we analyze the production of a secondary massive field \ lighter than the inflaton 
and we show that its production depends on the duration of inflation. In section XI we 
conclude and in the two appendices we relegate useful formulae for the adiabatic expansion 
and the dimensional regularization with cut-off. 



II. ANALYTIC APPROXIMATION 



We consider inflation driven by a classical minimally coupled massive scalar field. The 
action is: 



S = J d 4 xC = J d 4 x^f^g 



R 



-yd, 



ir vJ u (f) Til 



1x1 



(1) 



where £ is the lagrangian density and m is the mass of the field <fi. Further we consider the 
Robert son- Walker line element with flat spatial section: 

ds 2 = g^dx^dx" = -dt 2 + a 2 (t)dx 2 (2) 

The scalar field is separated in its homogeneous component and the fluctuations around 
it, <ft(t, x) = cf){t) +<p(t, x). During slow rollover the potential energy dominates and therefore 



K 



H z = — 



h m — 



,m 2 2 



(3) 



where k 2 = 8ttG = 87r/M 2 ]. The Hubble parameter evolves as: 



H = -y0 2 (4) 

On using the equation of motion for the scalar field: 

<j> + 3H<p + m 2 = (5) 

and neglecting the second derivative with respect to time we have 3H<p ~ — m 2 0, and from 
Eq. (W) one obtains: 



H ~ 

3 



(6) 



which leads to a linearly decreasing Hubble parameter and correspondingly to an evolution 
for the scale factor which is not exponentially linear in time, i.e.: 



H(t) ~H + H t a{t) ~ exp[H t + H Q -} 



(7) 



In Fig. [I] the comparison of the analytic approximation with the numerical evolution for the 
Hubble parameter H is displayed. 
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FIG. 1. We plot the numerical evolution of H (solid line) and the analytic approximation 
(dashed line). The initial condition corresponds to H{t = 0) ~ 9.2m and the time is in \ jm units. 



We now consider the equations of motion for the inflaton fluctuations <p in a rigid space- 
time (i.e. without metric perturbations): 



<f k + 3H k + u 2 k ip k = , 
where io\ = k 2 /a 2 + m 2 and the <£> k are the Fourier modes of the inflaton fluctuations, 



-ikx, „* 



(8) 



(9) 



As already emphasized in the introduction, exact solutions for scalar fields in an expand- 
ing universe are rare, and indeed we do not have exact solutions for Eq. (|j) with the time 
evolution given by Eq. (|7|). We therefore introduce an approximation scheme based on an 
analogy with de Sitter space-time, where exact solution for scalar fields with arbitrary mass 



and coupling to the curvature do exist. We introduce ifa 
dependence in -0 k as follows: 



,3/2 



The equation for ip^ is 



MC,h) with C = 4? 

an 



(fk and we split the time 



(10) 



^ 9 3 • 9 o 

— + m 2 - -H - -H 2 



0. 



(11) 



We now make the ansatz ^ k = Q^Zy^Q with /i, v and A functions of H . On expressing the 
first and second time derivatives as derivatives with respect to ((, H) and using H ~ 0, as 
follows from Eq. ffl), we obtain from Eq. flTT|) after a little algebra: 



8Z H 2 

ReSl C " +1 ~dt + ReS2 C+2 Zv + ReSa C " Zv + ° { JF : 



(12) 
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where we have neglected quadratic and higher order terms in H/H 2 . Indeed, in order to 
have a value of density perturbations compatible with observations, m is constrained to be 
0(1O~ 5 — 10 _6 M p i): from Eqs. (|3],|6]) one can see that working to first order in H/H 2 during 
slow-rollover (0 ~ fewM p i) is a good approximation. On considering H and ( as independent 
variables, the first term vanishes if Z u is a Bessel function of argument X( and index v. On 
requiring that the residual functions ReSj, i — 1 — 3 vanish individually, the parameters A, \i 
and v are determined to be: 

H H 

X = 1 ~Jp> ^=2H 2 (13) 
9 9 m 2 H 

»=-rw-*w < 14) 

Hence the general solution to Eq. (||) is: 

= ^ [AffW(AC) + BHf\\Q\ (15) 

where H^ 1 ' 2 ^ are the Hankel functions of first and second kind respectively, and A, B are 
time-independent coefficients to the order of our approximation. 

We note that in the de Sitter limit (H/H 2 0) the solution in Eq. fll5|) tends to the de 



Sitter solution ||10|jl9|]20| , since A = 1 , fj, — 0. 

On using Eq. the value for the index v in Eq. (|1^) corresponds to an exact scale 
invariant spectrum for the inflaton fluctuations (p, i.e. v = 3/2. We shall show numerically 
in section IV that this analytic approximation is very good for the relevant spectrum range. 
This numerical analysis agrees with a previous numerical estimate of the same spectral 



index [|21||. This scale invariant spectral index could seem a little surprising, since in de 



Sitter space-time, a mass term would lead to a spectrum, which is slightly blue shifted with 
respect to scale invariance (v < 3/2). To see this, it is sufficient to put H = in Eqs. 
(0,0). On considering H ^ (and negative), it appears to give a positive contribution 
to the mass term in Eq. (|TTD , instead it compensates the mass term in Eq. (0). The 
interpretation is the following: on considering a de Sitter stage in which H slowly decreases, 
a fluctuation freezes when it crosses the Hubble radius, with an amplitude determined by 
the value of the Hubble radius at the horizon crossing. However H decreases, therefore if 
kx > &2, this effect implies that the amplitude for the mode k\ is smaller than the one for 
the mode &2, since the latter crosses the Hubble radius first. This effect is a red tilt of the 
de Sitter scale invariant spectrum. For the case of slow-rollover in a chaotic inflationary 
model with a massive inflaton, these red and blue shifts exactly compensate, leading to a 
scale-invariant spectrum for the inflaton fluctuations ip in rigid space-time. 

III. QUANTIZED FLUCTUATIONS 

We now consider quantized fluctuations of the inflaton. This means that Eq. ([|) is 
promoted to an operator form: 
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ik x ; 



(16) 



where the bk are time-independent Heisenberg operators (also called time independent in- 
variants in |22 , 23| ). In order to have the usual commutation relations among the bk'- 

[6k, k>] = [bl bU = [6 k , bU = 6®(k - k') (17) 
one must normalize the solution to the equations of motion through the Wronskian condition: 

<^k <ft - = -3 ■ ( is ) 

This normalization condition yields to the following relation among the coefficients A, B of 
Eq. flTSp: 



\A\' 



\B\ 



71 

AH 



AC 2 "- 



(19) 



The fact that A, B depend on time should not surprise. In time- dependent perturbation 
theory, these coefficients, which would be time independent for exact solutions, acquire 
a time dependence just as the Wronskian of the solutions. In our case, this time 
dependence is consistent with the approximation, i.e. self-consistent to (including) order 
H/H 2 . 

The solution corresponding to the Bunch-Davies vacuum pO in de Sitter space-time, 
that is the adiabatic vacuum for k — > 00 during the slow-rollover phase, corresponds to 
choosing A = (irX/AH) 1 / 2 ^ 1 * , B = 0. With this choice, for X( — > 00, the solution (|15|) 
becomes: 



1 



y?k 



aV2k 



(20) 



Let us now discuss the behaviour for A£ <C 1. On using Eqs. ( JT~0|JT~3| ) one sees that \( <C 1 
implies k <C aH, i.e. wavelengths which are much larger than the Hubble radius. In this 
limit, the solution is PS ]: 



cp k ~ -1 



?{y) /ttA 



1/2 



Ti a 3 / 2 \4H 



(21) 



In order to compute expectation values of operators with respect to states in the time- 
independent invariant b p2|p3| basis it is useful to introduce the modulus of mode functions 

Xk 



l^k/v^l- The variable Xk satisfies the following Pinney equation: 

x k + 3Hx k + u\x k 



6,7,3 



a"x 



(22) 



We now rescale Xk = Pk/a 3 ^ 2 to eliminate the damping term and obtain: 



h + 



col - -H 2 - -H 



Pk 



Pk 



(23) 
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The general solution to Eq. (|23"D is given as a nonlinear combination of two independent 
solutions yi, i/2 to the linear part of Eq. (|23|) . From Eq. (|TT|) we can use the following 
solutions: 



y x = CMK) 

Therefore the solution to Eq. (p3|) is 

p k = (by* + My 2 2 + 2Ny lV2 



(24) 



(25) 



where the coefficients satisfy LM — N 2 = 1/W 2 , with W the (time-dependent) Wronskian 
of yi,y2- The choice of initial conditions for p k which corresponds to the adiabatic vacuum 
for k -> oo is N = and L = M = Att / \2H X 2 ^) . The solution for x k is: 



' f— V /2 AV 2 [^ (AC)+iV 2(AC)l 1/2 



a 3 / 2 V2i/; 

which coincides with the Bunch-Davies choice in the de Sitter limit 



(26) 



IV. NUMERICAL ANALYSIS 

In this section we present the numerical analysis of the time evolution of the ip modes. 
Besides checking the validity of the analytical approximation introduced in section II, this 
analysis is useful in order to understand how a natural infrared cut-off emerges in the 
problem, on assuming that inflation is not eternal in the past, but starts at some finite time. 
This infrared cut-off becomes relevant when the ip fluctuations are generated in an infrared 
state p7| , as occurs for v > 3/2 (see Eq. (|14D). 



First we want to analyze the properties of the spectrum of the inflaton fluctuations. We 
numerically evolve Eqs. (|3],[5D and Eq. (^2|). We present numerical data for an interval of 
comoving wavenumbers for which 1 < k/m < 10 5 at the initial time to (o(to) — !)• The 
initial conditions are 0(to) — 4.5M p i, 0(to) = for the inflaton. If we consider the vacuum 
state for each mode of the field fluctuations to be the initial condition, one has 

1 

Wo) - — - 1/2 

a3/2 (toW {t ) 

x k {t ) = (27) 

This fact can be easily seen in terms of the invariant operators introduced to quantize time 
dependent harmonic oscillators p2|-|2~I . 



As a second initial condition, we consider the limit for large k of the conditions in 
which correspond to setting the mass equal to zero. A third set of initial conditions is related 
to the adiabatic expansion in conformal time, (see later Eq. (|39|) ): 



a(t )nl /2 (t ) 

x k (t Q ) = -H(t )x k (t ) (28) 
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where Qk will be defined in Eq. (f40|). Let us note that for the last case the frequency becomes 
imaginary below a certain threshold and so we shall consider the region above it. In Fig. |2| 
we exhibit the three cases. The first two initial conditions lead to spectra practically equal 
for k of order m and above, the third set instead has a spectrum which joins the others at 
values of k of the order of Hq. 



o - 



-l - 



-2 - 



-3 - 




-10 12 3 

FIG. 2. We show on a logarithmic scale the r.m.s. of the infiaton fluctuations x^ obtained 

numerically for the three different sets of initial conditions. The spectrum is shown as function of 

Log(A;/m). 

The spectrum of the fluctuations, related to the initial conditions in (]27|) , are displayed, 
over a broader range, in Fig. [3[ 
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FIG. 3. We show on a logarithmic scale the r.m.s. of the infiaton fluctuations x^ obtained 
numerically (solid line) and and an asymptotic linear fit (dotted line). The spectrum is shown as 
function of Log(fc/m) and the linear fit is given by 0.7708 — 1.5 Log(fc/m), in agreement with [21]. 

Figs. and |3| display the spectrum at t = 10/m (for this case inflation lasts a period 
of time ~ 27/m). We note that this scale invariant spectrum extends only up to a certain 
scale, £, which is of the same order as the initial Hubble radius H(t ) (we have checked this 
by changing the initial conditions for the homogeneous mode of the infiaton). For comov- 
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ing modes smaller than t the spectrum oscillates and bends towards the initial conformal 
adiabatic vacuum state, as shown in Fig. (H). For all practical purposes we can therefore 
safely consider a scale invariant spectrum for k > I, where I = CH(to), with C a numerical 
coefficient 0(1). This numerical evidence favours the picture in which the amplification of 
the modes occurs mainly at the crossing of Hubble radius. Since m « H during inflation, 
this means that all the modes for which m <k % H(to), are not stretched by the geometry. 

Secondly, we wish to show how accurate the approximation fl2"6] ) is mode by mode. 
In order to do this we numerically solve Eq. (p3| ) and compare it to the approximation 
employed. The agreement is very good up to times very close to the end of inflation. 
In Fig. [| we show for the mode with k — 10 m the time evolution of the relative error 
— p^ pprox ^) / p^ um \ The larger k is the better the agreement. We are therefore 
allowed to use the approximation for the inflaton dynamics almost til the end of inflation. 
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FIG. 4. 

of 1/m). 



10 15 



20 



25 



Relative error (pjjT""^ — p" k 



(approx) )/p (num) ^ fc 



10 m) as function of time (in unit 



Let us now consider the correlator (p 2 ) (similar considerations are valid for other quan- 
tities bilinear in the field). We note from the spectral analysis that the integral over the 
modes can be split into two parts 



h 



(27T) 



d 3 k 



h 



2(2tt) 



d k Xu 



h 

471-2 



dk l\i oc i, ~\~ 



dk k 



2 2 



(29) 



Below the scale £ both the interval with the oscillations and the tilt around the scale 
invariant spectrum shown in Fig. 1 are included. The analytic treatment of the far infrared 
modes which contribute to the first integral would need an analytic approximation for the 
modes which at the beginning of inflation are outside the Hubble radius. This would amount 
to knowing the phase and the initial quantum states which precede the inflationary phase. 
Therefore we shall proceed by considering only the second integral and we shall neglect the 
first one in the far infrared. Even on neglecting the first integral which would require extra 



assumptions, the correct leading behaviour of the renormalized quantities is obtained [28 



We conclude this section by noting that other methods to deal with infrared states are 
present in the literature. Infrared states were studied by Ford and Parker |27]] for massless 
fields in Robertson- Walker space-times with a power-law expansion of a particular kind. By 
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matching an earlier static space-time with a space-time with a scale factor which expands 
in time with a power law, they noted that an infrared finite state cannot evolve to an 
infrared divergent state [|27|]. The same scheme was also used by Vilenkin and Ford |29j for 
the problem of massless minimally coupled scalar fields in de Sitter space-time. An earlier 
radiation dominated phase was matched to the de Sitter metric. Through this matching the 
infrared tail becomes suppressed leading to an infrared finite state. As we shall see both the 
calculations performed by eliminating the infrared tail (i.e. working with the cut-off) and 
suppressing the infrared tail (through the Bogoliubov coefficients obtained by the matching 
prescription) lead to the same result to leading order. Indeed, their physical motivation is 
the same: inflation is not eternal, but starts at a finite time. However, the two methods 
treat the infrared tail in a different way: the agreement to leading order implies that the 
relevant contribution to the correlator (y? 2 ) - and to the EMT - comes from intermediate 
modes, and not from the furthest infrared modes. 



V. THE ENERGY-MOMENTUM TENSOR 

The classical energy-momentum tensor (henceforth EMT) of inflaton fluctuations is: 



■-g al3 d a (pdf3(p - —<p 2 



(30) 



and its operator form is simply obtained by promoting tp to an operator as in Eq. 
When averaged over the vacuum state annihilated by b: 



&k|0) = 



(31) 



the energy-momentum tensor assumes a perfect fluid form because of the symmetries of the 
RW background [M: 



(Tfn,) = diag(e, a 2 p5, 



(32) 



where e and p are the energy density and the pressure density respectively. 

In the following we consider, according to the previous sections, v = 3/2, employing the 
dimensional regularization || to treat the UV behaviour. Therefore the integrands will be 
in 3 dimensions and the integration measure analitically continued in d dimension. 

The energy density is 



{Too) 



h 



2(2vr) d J\ k \>i 

h 



d d k 



4(27r) d J\k\>£ 



d d k 



k 2 

I . | o ' b i 1 2 2 i 1 2 

Wk\ + -z\<Pk\ +™> Wk\ 
or 



1 



6^2 



' k 2 



2 1 2 

+ m \x k 



(33) 



and the pressure density p, related to the space-space component of the EMT, is: 



P 



4(27r) d J\k\>e 



d d k 



4 + 



1 



a 6 xl 



+ 



'2-dk 2 
d a 2 



2 1 2 

m \x k 



(34) 
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On using Eqs. ( ^6|) and ( f73|) in Appendix B with a = d — 1, we may now compute the 
second part of the integral (p9f): 



2(27r)^| fc |>/ dk ^ 

-H 2 1 {2 -4 n2-2 — 

16tt 2 V 3# 2 /l V27T 1 



(35) 



d-3 



/2 




(36) 



Analogously, the energy and pressure densities in Eqs. (|33l , |3"4"|) , with the help of the 
formulae in Appendix B, are: 



e = (Too) 



16n : 



;H>{4 



m 



1 



H 2 V 27^/2 
2 



d-3 



nd—3 



r(i-d) 



m 



-l + 7 -21n2] + O (\)} + 0(d - 3) 



(37) 



P 



16n : 



m 



d-3 



m 



2 

l-7 + 21n2 



(38) 



The poles given by the negative values of the argument of the V function in Eqs. 
([36|,[37],[[8|) represent part of the ultraviolet infinities of field theory. 

VI. THE ADIABATIC SUBTRACTION 

In order to remove the divergent quantities which appear in the integrated quantities as 
poles in the V functions, we shall employ the method of adiabatic subtraction [[|. Such a 
method consists in replacing Xk with an expansion in powers of derivatives of the logarithm 
of the scale factor in Eqs. ([£||55|). This expansion coincides with the adiabatic expansion 
introduced by Lewis in |TjJ for a time dependent oscillator. 

Usually it is more convenient to formulate the adiabatic expansion by using the conformal 
time 7] 0] (drj = dt/a). We follow this procedure and write an expansion in derivatives with 
respect to the conformal time (denoted by ') for x/~- Then go back to the cosmic time and we 
insert the expansion in the expectation values we wish to compute. Adiabatic expansion in 
cosmic time and conformal time lead to equivalent results, because of the explicit covariance 
under time reparametrization [T2 . 



We rewrite Eq. ( p2]) in conformal time in the following way: 



(ax k )" + Vt\ (ax k ) 



(ax k f 



(39) 



where 
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k 2 + m 2 a 2 - -a 2 R 
6 



(40) 



and R is the Ricci curvature: 



R = 6- 



(41) 



The fourth order expansion for (<£> 2 ), the energy and pressure densities are therefore (as 
before v = 3/2, integrands in 3 space dimensions and the measure is analytically continued 
in d dimensions), using the expression in (|7Bp and the results of appendix A and B, 



h 



(4) 



16n' 



4 m 2 



T o 29 o 

— m 2 H H 

45 15 



am 
27V2 



d-3 



2 m 2 



r \2 2I + 9H 2 3 + 



(42) 



H(4) 



(Too) 



h 



(-1) 



16tT 2 " ^#2 V27T 1 /2 

119 m 2 
+ ~60 + IP 



d-3 



,d-3^d-3 



m d ^r(l - d) 



33 

h7 

10 ' 



(43) 



P(4) 



(r«> 



(4) 



167T 2 



# 4 {-4 



1 



H 2 V2vr 1 /2 



119 m 2 



d-3 



,d-3^d-3 



m d - d T(l - d) 



1309 



270 



-7 



+ 0(*O} + O(d-3) 



(44) 



VII. THE CONSERVED RENORMALIZED EMT 



On subtracting the adiabatic part given in Section [VT| from the bare integrated contri- 
bution given in Section [V| and taking the limit d — > 3 one obtains the finite renormalized 
expectation value for the correlator and for the energy-momentum tensor. 

The renormalized expectation value of (if 2 ) is therefore, neglecting terms of order 1/a 3 
and for a > H/m 



{(P 2 )ren 



("1) 



16n' 



i» 2 { I*- 




In a — In 



CH(t c 



7U 




41n2 - 



14 m 2 10 1 
'~9~JP + T ~~ ^ 2 



l m 2 + ^# 2 
45 15 



(45) 



Considering (|45), we note that for late times it resembles more a massless, than a massive, 
field in de Sitter space-time. This is a consequence of the scale invariant spectrum (|I~4| ) of 
inflaton fluctations (the same spectrum occurs for massless minimally coupled fields in de 
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Sitter space-time). The leading behaviour for (^ 2 )ren agrees for late times with the de 
Sitter result 0,^11 when H = and a(t) = a e HDst : 



2\DS 



REN 



4vT 2 



H DS t . 



(46) 



At earlier times, (^ 2 )ren is dominated by a contribution 0(H A /m 2 ), but then the time 
dependent piece 0(lna) takes over. We warn the reader about the massless limit taken at 
face value of Eq. (J4"5|). In the massless limit, as discussed in the appendix B, one has to 
use a different analytic continuation for the adiabatic part, related to the expression in (|77|) 
and the result for (^) 2 )ren will be finite, different from Eq. (^5|), but with the same leading 
contribution (|46|). However, this massless limit is just of academic interest, since for m = 
inflation would not happen. 

Even if both inflaton fluctuations for a massive inflaton and a massless minimally coupled 
scalar fields in de Sitter share the same scale invariant spectrum, the energy and pressure 
carried by fluctuations for these two cases are very different. For the latter case, a linear 
growth in time is present only for the correlator; the EMT does not contain the correlator, 
but only bilinear quantities less infrared than (p 2 (for a nice explanation of this difference 
see [0). For the case of inflaton fluctuations, the correlator appears directly in the EMT 
because of the nonvanishing mass. The kinetic and gradient terms should be smaller than 
the potential term, as for the massless minimally coupled case. Therefore for the case of 
inflation driven by a massive inflaton, the EMT of inflaton fluctuations should grow in time. 
This is what we shall show in the following. 

On subtracting Eq. ([43]) from Eq. fl3~7|), the renormalized energy density €ren is: 



CREN 



(Too) REN 

h 



167T 



(Too) — (Too) (4) 
CH(t ) 



I h 2 



In 



711 



lna(t) 



119 
~60~ 



2 

m 



— -21n2 
10 



O 



Similarly, by subtracting Eq. ( fill) from Eq. ([38)), the renormalized pressure density Pren is: 

(^ii)REN _ (Tu) — (Tu)(4) 

119 



Pren 



I 7T2 



16tt 2 



ln^^-lna(t) 



m 



Jl 2 



1039 
~270~ 



+ 21n2]+o(l)} (48) 



We note that this result does not agree for H = (and therefore for m = because of 
Eq. (H)) and £ = with the de Sitter result |]5|j30]| obtained with the Bunch-Davies vacuum: 



thDSBD 
1 fiv REN 



9jw 
64vr 2 



m 



m 2 + (f 



6' 



)R 



,.3 . .,3 N . 12m 2 
^2 +")+^2 _Z/) " 



K ^ 6 J 18 2 V ^ 6 J 2160 



(49) 



where R is the curvature in de Sitter (R = 12ifp S with H^s as the Hubble parameter in de 
Sitter space-time) . Indeed, the limit of Eq. (|49[) for vanishing mass and coupling is 

finite and is [fEll : 
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thDSBD 
/ii/ REN 



61^s 
960tt 2 



Instead, the result (f4"T|,4*B) agrees with the Allen- Folacci [13] result for m = 0: 



t-iDS AF 
J /i^ REN 



960vr 2 



9nv 



(50) 



(51) 



However, the renormalized EMT of inflaton fluctuations in Eqs. (H^H) grows in time, 
as the logarithm of the scale factor. This feature is due both to the fact that the inflaton is 
massive and to the infrared state in which its fluctuations are generated. For a test field in 
de Sitter space-time a linear growth in time of the EMT is possible, only for m 2 + £R = 0, 
with m and £ both different from zero. 

The renormalized EMT in de Sitter space-time in the Bunch-Davies vacuum 

(H) and 

in the Allen-Folacci vacuum corresponds to a perfect fluid with an equation of state w = 
p/e = —1, which is identical to the background driven by a cosmological constant. The 
conservation of the renormalized EMT (f49|,|51~D is direct consequence of its symmetries: 



^o/ S REN fiV 



v 1 nv REN 







(52) 



"9^ = 0. 

The renormalized EMT in an inflationary stage driven by a massive inflaton corresponds 
to a perfect fluid, but with an equation of state which differs from —1 by terms 0(m 2 /H 2 ), 
as one can see from Eqs. ( |47l , [48| ) . 

The derivation of the conservation of the renormalized EMT in chaotic inflation is also 
straightforward. The renormalized EMT is conserved consistently with the approximation 
used, i.e. to the order 0(m 2 /H 2 ). The conservation can be easily checked mode by mode, 
i.e. by considering the covariant derivative inside the integrals in k in the difference between 
the bare value and the fourth order adiabatic value and using the equations of motion for the 
field modes (§). The conservation of the final renormalized value EMT can also be checked 
by inserting the expressions given by Eqs. (|47|) and ( 



m 



de 



REN 



dt 



+ 3H (e RE N + Pren) = . 



(53) 



and retaining only the terms up to and including 0(m 2 /H 2 ). 



VIII. BACK-REACTION ON THE GEOMETRY 

We now discuss the back-reaction of the amplified ip fluctuations on the geometry. 

We consider the back-reaction equations perturbatively: we evaluate the higher order 
geometrical terms |J generated by renormalization as their background value and we do 
not use them to generate higher order differential equations. We note that, in accord with 
the approximations used, many higher order derivative terms are implicitly absent in Eqs. 
( |4?| , fl8|) , since they would be of higher order in powers of H/H 2 and because H ~ 0. We 
estimate the back-reaction effects without changing the structure of the left hand side of 
Einstein equations. Hence the back-reaction equations we consider are: 



14 



H 2 
H 



8nG \(j) 2 m 2 2 

— + —<p + e REN 

-4:71 G((}) 2 + 6 REN + PREN) ■ 



(54) 
(55) 



The main point is that the energy and pressure of inflaton fluctuations grows in time 
as the logarithm of the scale factor, while the Hubble parameter driven by the background 
energy density decreases linearly in time. However, the approximation we use, i.e. Eqs. 
(§,[7|), is valid for a certain time interval, At given by: 



Ho 

2 



At ~ ^ (56) 
Therefore the term which grows in the renormalized EMT will saturate at the value: 



e REN (At) ~ -p REN (At) 



in 



rr2 tj2 
2 H H 



(57) 



The term ( (5T|) is larger than the Allen-Folacci value ~ if 4 in Eq. ( |5TD and has the 
opposite sign. Therefore, the energy density of fluctuations starts with a negative value and 
changes sign to a positive value when the logarithm takes over. This behaviour is shown in 
Figs. [| and ||. 

Such a value leads to a contribution to the Einstein equations of the order H 2 Hq/M 2 v 
The importance of back-reaction is therefore related to the ratio Hq/M 2 v If inflation starts 
at a Planckian energy density then back-reaction during inflation cannot be neglected. 

A sligthly different conclusion can be reached on considering the variation in time of 
the Hubble parameter, i.e. Eq. (|55|) . The important point to note is that the leading 
contribution in 6ren,Pren, i-e. the terms ~ H A and ~ m 2 H 2 log a, do not contribute to 
H since these terms have an equation of state PrenAren = — 1- The contribution of the 
inflaton fluctuations to H is therefore of order m 2 H 2 /M 2 1 , which is suppressed with respect 
to the classical value by the factor H 2 jM 2 x . 
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5. Time evolution of the renormalized energy, where the cosmological time is in units of 
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FIG. 6. Time evolution of the magnitude of the back-reaction 8-7T G €ren /(3H 2 \ eREN= o), where 
the cosmological time is in unit of 1/m. 



IX. COMPARISON WITH THE SLOW-ROLLOVER CALCULATION 



We now compare the approximation which lead us to Eq. ( [L5D with the slow-rollover 
technique introduced by Stewart and Lyth QTSJ . The latter technique was developed directly 
for scalar and tensor metric perturbations (jTj| , and not for field perturbations in rigid space- 
times, as treated here. However, the equation for gravitational waves differ from Eq. (^|) 
only by the presence of the mass term m. Therefore, as a first check we note from Eq. fll4"|) 
that 



GW 



9 
4 



H 
''W 2 



3 
2 



H 



(58) 



where the second relation holds when H/H 2 is small. The value i^gw i n Eq. 

(41) of 



5q ) coincides 

On the other hand, if one applies the Stewart-Lyth 



with the index \x of Eq 
procedure to Eq. (H) one would obtain an index for Hankel functions which is precisely v in 
Eq. fllj. 

A natural question is to ask whether inflaton fluctuations are generated in infrared states 
also for other models of chaotic inflation. Analytic approximations, such as the one presented 
in Section 2, are very difficult to obtain. However, since our calculation agrees with the slow- 



rollover result [|T8| for inflaton fluctuations in rigid space-times and gravitational waves in 
the case of a massive inflaton, one can use the latter technique to estimate the spectral 
index of fluctuations in a generic chaotic model with potential V(<f>) = X(p n /n (here A has 
the dimensions of a mass elevated to the power 4 — n). We follow Stewart and Lyth and we 
study the equation: 



(aip k )" + [k 2 + ml s a 2 {a<p v 



0. 



(59) 



where m 2 s = V,^ = X(n — l)0 n 2 . On using a{rj) = —l/[Hr](l — e)), after some algebra one 
obtains the following spectral index for the constant mode: 
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1 l-2e + 2e/n 
v = — I 

2 1 -e 

where the slow-rollover parameter e is taken as constant and is defined as: 



(60) 



H 



n 



e = — 



H 2 167TG0 2 



(61) 



The slow-rollover approximation is better for large values of the inflaton. The result ( |60| ) 
agrees for n = 2 with the quadratic case. However, for n > 2 the spectral index v is smaller 
than the critical value 3/2, leading to inflaton fluctuations which are infrared finite. 

Although the approximation presented here gives the same result as the slow-rollover 
technique to the lowest order |18| for gravitational waves and inflaton fluctuations in rigid 
space-times for a massive chaotic model, it will not be so when one includes metric pertur- 
bations WK . 



X. MODULI PRODUCTION 

Let us now discuss the quantum production of a light scalar field x i n this model of 
massive chaotic inflation. By the term light, we mean a scalar field with a mass M, which 
is smaller than (or equal to) the inflaton one, m. We assume a vanishing homogeneous 
component for x- Therefore, for x the approximation of a rigid space-time is correct. 

With M < m, we see from Eqs. (||,14) that the index v x for the Hankel functions, which 



are involved in the solutions for the x modes, is larger than 3/2. The leading contribution 
for the renormalized value of (x 2 ) is: 

(xVen „ m 

where (3 is a numerical coefficient. This formula generalizes Eq. ([4q ) to the case of a mass 
M < m. The case M = m represents a limiting case, for which a logarithm appears. 
Analogously, in de Sitter space-time, a limiting case is also present for v = 3/2 Hl2"| , and a 
logarithm of the scale factor appears instead of a power. 
For a massive % the following relation holds: 

3 



where z/qw is defined in Eq. (E3|) and, owing to the smallness of the ratio m 2 /H 2 , v x is 
very close to 3/2. However, the growth in time of (x 2 )ren is more rapid than the growth of 

This fact can be checked numerically: for example, on taking the initial condition for 
inflation already used in the previous sections for the other numerical checks, one finds for 
M = 0, v x ~ 1.5043 and for M = 0.5m, v x ~ 1.5032. Such a value is time independent 
after a very short transient phase needed for the modes to freeze and corresponds to the 



one given in (|58 ) for H computed for a time very close to the beginning of the inflation. 



We can therefore substitute H{t ~ 0) ~ H$ in place of H(t) in v. We have checked this 
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relation for different initial Hq. This fact is in agreement with the fact that the approximate 
solution, which we employed for the inflaton field, can be used for other fields with different 
masses, but only at the beginning of the inflation period. Thus we can control the spectral 
behavoiur of the moduli field, but not its normalization. 

To study the behaviour in (|62| ) it is convenient to rewrite the evolution of the scale factor 
a and its exponent (2^ x — 3) as 



a(t) 



exp - 



3H 2 



H 2 



2m 2 -M 2 

2z/ v — 3 ~ k 

x 3 H 2 



(64) 



Therefore the main result (let us only write the dominant contribution) is that for M < m 
and at the end of inflation 



(x 2 )ren ~ 



H 2 H 2 



exp 



H 2 ~H 2 



1-Ml 



1 



M 2 



'2\ . H 2 



REN 



HI 

2 



m 



(65) 



and for the renormalized EMT associated with the \ field 



e REN ~ 



——H 2 H 2 

167T 2 111 2 " 



exp 



H 2 



M 2 



M 2 
m 2 



(66) 



The results in Eqs. (57|,65|, |66f ) are very interesting. They show that the production of a 
scalar field x with mass M smaller or equal to the inflaton mass m depends on the duration 
of inflation and is larger than the usual extrapolation of the de Sitter result {(x 2 ) ~ H 4 /m 2 
and e x ~ (M 2 /m 2 )H A ). This result implies that the quantum production of light scalar 
fields in chaotic inflation with a mass term is even greater than expected on extrapolating 
the de Sitter result, and depends on the duration of inflation, as is also stated in pl| , |35| . 
Again, if inflation starts at a Planckian energy density, the back-reaction of light scalar fields 
cannot be neglected during inflation. 

We note that the factor in curly brackets in (|65| , |66|) is larger than the logarithmic term 
one has for the v — 3/2 case, by a factor of 2. 

We see that on computing the exact numerical solution of the Pinney equation for the 
modes of the inflaton with mass m and for a moduli field with M = 0.5m, one obtains for 
the same momenta amplitudes more than one order of magnitude larger for the latter at the 
end of inflation, again in agreement with [ 21 1 . We therefore expect an enhancement of 2-3 



orders of magnitude for its backreaction with respect to the inflaton case. 



XI. DISCUSSION AND CONCLUSIONS 

We have computed the renormalized conserved EMT of the inflaton fluctuations x) 
in rigid space-times during the inflationary stage driven by a mass term. The method of 
dimensional regularization has been applied by using an analytic approximation valid during 
the slow-rollover regime. All the results agree with the Allen-Folacci results for of a test 



field in de Sitter space-time [14], in the limit for which the Hubble parameter is constant 



(which is also the massless limit because of Eq. 
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We find that the EMT of inflaton fluctuations grows in time. The reason for this be- 
haviour is that chaotic inflation driven by a massive scalar field produces a scale invariant 
spectrum of fluctuations even if the field is massive. This effect is due to the decrease of the 
Hubble parameter during the slow-rollover regime. 

In de Sitter space-time, the renormalized EMT of a quantum field grows linearly in 
time only if m 2 + £-Rds = with m and £ different from zero [12]. A massless minimally 
coupled scalar field in de Sitter space-time, characterized by a scale invariant spectrum of 
fluctuations, leads to a correlator which grows in time. However, only bilinear quantities 
less infrared than the correlator appear in the EMT, and therefore the expectation value of 
the EMT of a massless minimally coupled scalar field is constant in time ||12| . In massive 
chaotic inflation, inflaton fluctuations are generated with a scale invariant spectrum. Since 
the correlator appears directly in the EMT because of the nonvanishing mass, then the 
renormalized EMT grows in time just as the correlator does. 

We find that the growth of the EMT of inflaton fluctuations during slow-rollover leads 
to a positive energy density which reaches a maximum value 0(H 2 Hq), where H is the 
Hubble radius at the beginning of inflation. This value exceeds the usual value 0(H 4 ), 
which is of the same order of magnitude as the conformal anomaly. These values also show 
that back-reaction effects cannot be neglected if inflation starts at Planckian energies, i.e. 
at H ~ M p \. If inflation started at Planckian energies, although the contribution of the 
terms 0(H 2 Hq) and of the conformal anomaly would be of the same order of magnitude, 
we think that the two contributions could be different because of the different signs and of 
the different behaviours in time. 

In this model of chaotic inflation, we have also analyzed the geometric production of an 
additional field x with mass M smaller than the inflaton mass m. Of course, on considering 
the normalization, we have found that e x > e v ~ H 2 Hq. This result implies that the 
quantum production of light fields depends on the duration of inflation and it is greater 
than expected on extrapolating the de Sitter result (in de Sitter e x ~ H 4 ). As in the case 
of inflaton fluctuations, the energy density of light scalar fields could be comparable to the 
background one at the end of inflation, if inflation started at Planckian energy densities. 
Also in this case, the back-reaction of \ fluctuations does not appear to be negligible during 
inflation. 

One may then ask whether this behaviour of the back-reaction due to the fluctuations in 
rigid space-time is common to other inflationary chaotic models. Analytic approximations, 
such as the one presented in Section 2, are very difficult to obtain. However, since our 
calculation agrees with the slow-rollover result |18|] for the massive case, we have used the 
latter technique to estimate the spectrum of inflaton fluctuations in rigid space-time for 
a generic inflaton potential V(4>) = \<ft n /n. We have found that for n > 2 the inflaton 
fluctuations are generated in an infrared finite state, leading to a back-reation which does 
not increase in time. However, we think that we must address the problem while including 
metric perturbations, in order to fully understand this issue. It is known that chaotic 
inflationary models predict a spectrum of curvature perturbations which is red tilted JTB| - 
i.e. with a spectrum more infrared than the scale invariant one -, a result which does not 
hold for field perturbations in rigid space-time, as we have shown. Since infrared states 
could lead to a correlator which grows in time, the possibility exists that a back-reaction 
growing in time is common to all the chaotic inflationary models once metric perturbations 
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are included. 

Other important issues are whether an eventual self-consistent scheme to include the 
back-reaction would prevent the development of infrared states. The effect of the self- 
consistent inclusion of back-reaction effects on the spectrum of fluctuations during inflation 
is, to our knowledge, an issue still to be fully explored. It would be interesting also to 
investigate the effect of the inclusion of the higher order terms |3f| in the back-reaction 
equations ( p4] , po] ). Obviously, the Starobinsky model is a surprising example of the 
importance of higher order terms. 
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XII. APPENDIX A: THE ADIABATIC FOURTH ORDER EXPANSION 



From Eqs. (p9|,^0|,|41 ) one obtains the expansion for Xk up to the fourth adiabatic order: 
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(67) 



where Qk is defined in Eq. 



and €2 , 64 are given by: 

in: 
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(68) 



The solution in Eq. fl6"7|) must be expanded again since the Ricci curvature is of adiabatic 



order 2. Therefore is: 
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(69) 



where c = a 2 and 
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S fc = {k 2 + a 2 m 2 ) 1 / 2 
1 Si 3 T! 2 



e 2* — TTTvT + 7 



15 

4S 



3 6 2 



4S| 



1 1 

4S 



2 6 2 



(70) 



XIII. APPENDIX B: DIMENSIONAL REGULARIZATION WITH CUT-OFF 



We start with if 2 as an example. According to the discussion following Eq. (|29|), we 
neglect the infrared piece of the integral since it gives a small finite part. Therefore the 
relevant integral 



<^ 2 > 



1 h ?7r 3 / 2 r+°° 
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ft A 7T 27T 3 / 2 /■+ 
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(71) 



on extending it to (i-dimensions (integrands in 3 space dimensions and analytic continuation 
of the measure to d dimensions) : 



<^ 2 > 



h A 7r 2ir d / 2 r+ 



dk k 
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J2 (xk\ xk; 

v \aH u \aH. 



(2n) d a 3 4HT(d/2) 

The two-point function can be computed by using the following integral [37]: 
dxx a [J 2 (x) + N 2 (x)} = 



(72) 
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(73) 



With u — 3/2 and a = d — 1, and on using: 
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the expectation value for (<£> 2 ) is: 
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(75) 



Similarly, the integral used for the fourth order adiabatic quantities in <i-dimensions is 
(integrands in 3 space dimensions and analytic continuation of the measure to d dimensions): 
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Let us note that on taking the massless limit one can analytically continue the hypergeo- 
metric function and after straightforward calculations one gets 
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(77) 



which is of course the result one would obtain setting m = from the beginning. Therefore 
all the massless singularities in ( fT6| ) correctly cancel. We note that in the massless limit the 
analytic continuation misses the UV divergencies stronger than the logarithmic ones. 

On again considering the case m ^ such that a(t) > H/m and using the result ([TT|). 
which also holds also for 2-F1, one obtains, to the fourth adiabatic order, for (<£> 2 ): 
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